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Abstract 

This paper is devoted to the analysis of an abstract formula describing quantum adiabatic 
charge pumping in a general context. We consider closed systems characterized by a slowly 
varying time-dependent Hamiltonian depending on an external parameter a. The current oper- 
ator, defined as the derivative of the Hamiltonian with respect to a, once integrated over some 
time interval, gives rise to a charge pumped through the system over that time span. We de- 
termine the first two leading terms in the adiabatic parameter of this pumped charge under the 
usual gap hypothesis. In particular, in case the Hamiltonian is time periodic and has discrete 
non-degenerate spectrum, the charge pumped over a period is given to leading order by the 
derivative with respect to a of the corresponding dynamical and geometric phases. 



1 Introduction 

1.1 Motivation 

Many physical systems of interest can be described by means of a time dependent Hamiltonian 
reflecting the action of external agents on the system or taking into account the variations 
of its environment, in an effective manner. In such generality, little can be said about the 
evolution of the system. However, when the Hamiltonian is a slowly varying function of time 
with respect to some fixed relevant time scale of the system, the adiabatic theorem of Quantum 
Mechanics provides a very useful tool to describe the evolution in an approximate way under 
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certain hypotheses, see e.g [S][32]. The mathematical circumstances under which an adiabatic 
theorem of Quantum Mechanics can be proven are diverse, starting with the well known gap 
hypothesis in the spectrum of the Hamiltonian, see e.g. [191 1231 [4], which will be enough for 
our purposes. Note, however, that higher order approximations 24, 25, 18 , generalizations to 
situations where the gap assumption fails [9j[2] or where self-adjointness does not hold [281 117] 
or to a space-time setting [33] have been carried out over the years; see also the review [15] and 
references therein. 

In several physical systems, the Hamiltonian also depends on some external parameter, 
a £ R, associated with an instantaneous charge current. For instance, this is the case for 
models used in the study of the Quantum Hall effect where the electric current is monitored 
by a magnetic flux, a, through the sample which appears as a parameter in the Hamiltonian 
H a (s). The instantaneous current operator is given by d a H a (s) see e.g. [261 14], A similar 
phenomenon holds in models of mesoscopic physics where the current across the device they 
describe is driven by a phase difference. Examples are superconducting systems used as Cooper 
pair pumps [31] Q] [T3] [£1 [29] [30] \W\ . 

Accordingly, for an initial state ip, the charge pumped through the system in that state 
between time and r is given by 

<Q Q (r)) v , := / (^7 a (T')VI(a Q // Q (r , ))^7 a (r')V>^ , , 
Jo 

where U a (r) is the evolution operator. It is often true in applications that the time variation 
of the hamiltonian is periodic, in which case one considers the charge transported over a time 
period. For a slowly varying Hamiltonian H a with an isolated part of spectrum and an initial 
states tp in the corresponding spectral subspace, the pumped charged {Qa(i~))^ can be computed 
by making use of an adiabatic approximation. This yields the starting point of the analysis of 
the geometrical and topological properties of this quantity. In particular, when the dependence 
of the hamiltonian is periodic in both the time and parameter a, it is well known that the 
pumped charge over a cycle suitably averaged over a is quantized and related to topological 
indices associated with the spectral projector, see e.g. [26[[5lf4]. However, in certain mesoscopic 
devices, the observed charge transfers do correspond to {Qc«(t))^ for fixed values of a [10] . 
which is also known do differ significantly from the average value over a, [3]. 

Note that in contrast to other mathematical studies, see e.g. [3], the time-dependent adia- 
batic pumps we consider are not open systems in the following sense: their dynamics for frozen 
times is not characterized by an explicit scattering matrix between infinite reservoirs. Therefore, 
no formula is available to determine the charge transport in terms of matrix elements of the 
instantaneous scattering matrix [111 112] . 

The purpose of this paper is to provide a mathematical derivation of the adiabatic pumped 
charge through the system in an abstract setting that covers the physical situations described 
above. 

We focus on the derivation of a controlled approximate expression for (Q a (r))^,, the charge 
pumped through the system over a time span [0, t] for fixed values of a up to vanishing correc- 
tions in the adiabatic regime. We work in a quite general framework, with arbitrary dependence 
of the Hamiltonian on time and a, assuming only the existence of an isolated part of spectrum 
in H a (r). This allows to recover as particular cases the time-periodic situations and, for an 
isolated eigenvalue in the spectrum, expressions for the charge pumped over a period used in 
several contexts, involving the geometric phase or its non-abelian version [10] . 

Moreover, we hope the present analysis can provide an introduction to the mathematical 
methods used in the analysis of adiabatic phenomena beyond the matrix case. 
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1.2 Setup and Results 



Let us describe more precisely the mathematical setup and results we prove in the next Section. 
The adiabatic regime is characterized by a dimensionless time scale 1/e, which is assumed to 
be long, i.e. e <g 1. The Hamiltonian is assumed to depend on e as follows: H a — _ff a (er). 
This makes the Hamiltonian slowly varying since it takes a time of order 1/e for it to change 
significantly. Introducing the rescaled time variable t — er, of order one, the corresponding 
evolution operator U a also depends on the small parameter e and satisfies 



Accordingly, for any initial state ip, the charge pumped over a time interval [0, r = t/e] now 
reads 



which is of order 1/e. From here on, s is consistently the rescaled dummy integration time 
variable. In order to capture the leading term, up to a vanishing error in the adiabatic limit 
e — > 0, it is necessary to compute the next to leading order approximation of the evolution 
operator in e. We do this in Theorem l2.2l under the sole gap hypothesis, after having recalled the 
usual leading order adiabatic approximation in Theorem l2.ll Then we focus on initial conditions 
that belong to the spectral projector corresponding to the isolated part of spectrum by the gap. 
The corresponding general expression for the charge operator is provided in Proposition 12. II A 
far more tractable expression is obtained for an isolated eigenvalue, of arbitrary degeneracy, as 
Corollarv l2.ll In case the Hamiltonian is time periodic oft period one, we recover in Corollary 1 2. 2 1 
the expressions used in |10] , The simplest instance being for a simple isolated eigenvalue E a (t), 
such that H a (t)^(t) — E a (t)ip(t). The charge pumped over a cycle for the initial condition 
= ip{0) then reads 



Here the first term is the usual dynamic contribution, whereas the second one is of geometric 
nature [32U34|, j3 a being the familiar geometric phase [51 [7]. In case the eigenvalue is degenerate, 
the matrix valued non-abelian generalization of this quantity, -B(l), to be used is of course also 
described in Corollary 12.11 

The last Section of the paper is devoted to the study of a family of examples used in the 
physical application described in [10] . The Hamiltonians considered display a permanently 
degenerate isolated eigenvalue for which explicit computations of B(l) can be performed. In 
case the dimension of the corresponding degenerate eigenspace is two, which corresponds to the 
applications considered in [TO], we also provide a geometric interpretation of B(l). 



In this section, we provide the abstract rigorous mathematical analysis behind the formulas 
used in study of the physical phenomenon of adiabatic charge pumping. We feel such a rigorous 
analysis is useful because of the presence of a variety of formulations of geometric adiabatic 
charge pumping in the litterature which have similar features. This allows to make clear under 
which hypotheses we work. Also, we believe the analysis is interesting in itself because it applies 
under very general conditions and it might be of use in different frameworks. 

Let T-L be a separable Hilbert space and H a (t) be a bounded self-adjoint operator on H. In 
order not to obscure the analysis by side issues regarding the technical difficulties related to the 



ied t U a {t) = H a (t)U a (t), U a (0) = I. 
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use of unbounded operators, we stick to the bounded case. For the same reason, we assume the 
parameter a is real valued. We work under the 

Regularity assumption: The map (t, a) h-> H a (t) is C 3 in the norm sense, as a bounded 
operator valued function, with (t, a) £ [0, 1] x [0, 1]. 

Let U a (t) be the solution to 

ieU a (t) = H a (t)U a (t), U a (0)=I, (1) 

where we rescaled time for convenience. Since C(W) is a Banach space, it follows from the 
general theory of differential equations, see e.g. |21j Ch. VI, that the solution to (fTJ is as 
regular in (t, a) as the Hamiltonian is, i.e. 

(i,«)^M')£C s ([0,l] x [0,1], C{H)). (2) 

The variable a is a parameter whose variations monitor fluxes or currents in the time de- 
pendent physical device described by H a (t). The current operator being defined as d a H a (t), 
if ip £ H is a normalized initial state the average charge pumped by the system between the 
physical times and t/e is equal to 

(Q a (t))i = - [*{U a (a)il>\(d a H a (a))U a (s)il>)d8. (3) 
e Jo 

This leads us to study the charge operator Qa(t), whose matrix elements carry the physical 
interpretation, defined by 

Q a (t) = - f U-^s^H^U^ds. (4) 
e Jo 

The geometrical properties of another average charge operator in the adiabatic limit have been 
investigated in [5] under the assumption that (t, a) >-> H a (t) is periodic in both variables. The 
average there is taken both over the period of the time dependent Hamiltonian and over the 
flux variable a. By contrast, we analyze the adiabatic behaviour of the charge operator for fixed 
values of a and t. In case of a time periodic Hamiltonian, the geometric content of the charge 
operator over a period is elucidated. 

As already observed in [4] for example, an alternative exact expression for Q a (t) reads as 
follows: 

Q a (t)=iU- 1 (t)d a U a (t). (5) 
It is a consequence of the regularity of U a (t) and the computation 

iedtiU-itydcUait)) = -U-\t)H a {t)d a U a {t) 

= U~\t){d a H a {t))U a {t), (6) 
with d a U a (0) = 0. Note that we also deduce from ([5]) the general estimate 

d a U a (t) = 0(t/e), (7) 

which is sharp in the scalar case. 

We proceed by considering the adiabatic limit e->0 under the familar 
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Gap hypothesis: Assume the spectrum a(H a (t)) consists in two disjoint parts a(H a (t)) — 
a a (i) U o^r(i) such that 

inf dist (cr„ (*) , ovT(t) ) = g > 0. (8) 

(t,a)S[0,l] 2 

Let P a (t) be the spectral projector of H a (t) associated with a a (t) by means of the Riesz 
formula 

Pa(t) = -^-.J(H a (t)-z)- 1 dz, (9) 

where 7 is a loop in the complex plane encircling a a (t) only, which is locally independent of 
(t, a) and let P a (t) = I — P a (t) be its complement. These projectors are as regular as the 
Halmiltonian is and moreover satisfy for any (t, a) 

P a (t)P a (t)P a (t) = P«(t)(d a Pa(t))P a (t) = 0, (10) 

as easily seen by differentiating the identity P a (t) 2 = P a (t). 

We introduce two unitary operators whose product will approach the true evolution in the 
adiabatic limit. Let Wait) and <& a (t) be defined by 

iW a (t) = K a (t)W a (t), W a (0)=I, (11) 
«■*„(«) = W~\t)H a (t)W a (t)$> a {t), $ a (0) = I, (12) 

where 

K a (t) = i[P a (t),P a (t)}. (13) 
It is a classical fact that the following relations hold, 

W a (t)P a (Q) = P a {t)W a (t) (14) 

$ a (t)Pa(0) = P„(0)*a(t). (15) 

The first identity is proven by showing that both operators satisfy the same differential equation 
with same initial condition, exploiting the relation 

P ol {t)K a {t)P a {t) = Q, (16) 

which is a consequence of (|10|) . see e.g. [T5]. The second identity follows from the fact that, by 
construction, the generator of & a (t) commutes with P Q (0). 

We have the 

Theorem 2.1 (Adiabatic Theorem) Assuming the Regularity and Gap hypotheses above, we 
have for any (t,ot) £ [0, 1] x [0, 1], 

U a (t) = Wc,(t)$ a {t) + 0(e), (17) 

where the error term is uniform in a £ [0, 1] . 
Remarks: 

i) As stated, the Theorem dates back to [23], [24] and was generalized by [4] to the unbounded 
case. The fact that the error term is uniform in a £ [0, 1] is a straighforward consequence of 
our Regularity Assumtions and of the Gap hypothesis which is uniform in a £ [0, 1] . 

ii) In case a a (t) consists of a single eigenvalue E a (t), the Theorem says 

U a (t)P a (0) =e-i^ E ^ a)ds P a {t)W a {t) + 0{e). (18) 
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This is the statement proven by Kato in |19| . 

ii) Further assuming E a (t) is nondegenerate and the time-dependent Hamiltonian is periodic in 
time, of period 1, we get the geometric or Berry phase out of this formula as follows: 

Let <p Q (0) = -P Q (0) ip a (0), be a normalized eigenvector associated with E a (0). Then, property 
(Tl4)l implies that V(a,i) G [0, l] 2 

<p*(t) = W a (t)<p a (0) satisfies H a (t)cp a (i) = E a (t)<p a (t). (19) 

and (TTTj) together with (TI5J) yield 

(<p a {t)\<p a (t)} = 0. (20) 
By periodicity, ip a (0) and faiX) differ by a phase only 

= e~ l % a (0), (21) 

where f3 a is the geometric or Berry phase. Therefore we finally get 

!/«(]>« (°) = e ~i^ BQ(3)ds e- l/3 > a (0) + O(e). (22) 

This shows concretely that the operator W a it) carries the geometrical content of the adiabatic 
Theorem. Note finally that if 1 ij) a (t) is another choice of instantaneous normalized eigenvec- 
tor associated with E a (t) which is periodic in time, these vectors differ from tp a (i) by a phase 
for any time 

Va (t) = e-^ [t) ^ a (t). (23) 
Using (|20|) . we get an explicit expression for b 

b a (t) = M0) ~i [ (Ms)\Ms))ds. (24) 
Jo 

Therefore tp a (l) = e~ $ ^ a (s)l * a (s)>d Vc<(0) and 

P*=-i [ {il> a (8)\i> a (s))ds. (25) 



Coming back to the charge operator, we see that in order to estimate Q a (t) up to errors of 
order O(e) vanishing in the adiabatic limit, we need to control the evolution to order 0(e 2 ), see 
©. This can be achieved as follows, see e.g. [24], [4], [18], [25], 

Let 

Hi 1 \t) = H a (t)-eK a (t), (26) 
which satisfies the gap assumption ([8j) for e small enough since K a (t) is bounded i.e. 

a(Hi 1) (t)) = ai 1) (t)U^P(t). (27) 

Hence we can define the corresponding spectral projectors P« (t) by ([9| and (t) = I— (t)- 
By perturbation theory we have for e small enough 

{H<£\t) - z)" 1 = (H a (t) - z)- 1 + (H a (t) - z)-\K a {t){H a {t) - z)- 1 +0{e 2 ) (28) 

where the remainder term is C 2 in (t, a) and of order e 2 , as a uniformly convergent Neuman 
series. The same is true for the perturbed projector expressed using ([9]) with the same path 7 
for e small enough, 

P£\t) = P a (t) - eK a (K a )(t) + 0(e 2 ), (29) 



G 



where, for any bounded operator B, 

n a {B){t)= ^- j {H a (t) - z^BiHait) - z)- x dz (30) 

with 7 is a loop encircling a2~\t), which can be chosen locally independently of (t, a). Let us 
also note here for future reference that 

P a (t)TZ a (K a )(t)P a (t) = 0, (31) 

as a consequence of the fact that the resolvent and the spectral projectors of H a (t) commute 
and of {T6]). We set 

K^(t) = i[Pi 1) W,P ( < 1) W] = K a (t) + 0(e), (32) 

where, again, the error term can be differentiated without harm. We define Wa (t) and (t) 
by 

iW^(t) = KM(t)wW(t), Wi 1) (0)=I, (33) 

ie^(t) = w£ ) ~\t)( K H£\t) + eVl{K a {t))w£\t)^\t), (34) 

*L 1) (o) = i, 

where, for any bounded operator B, T)\ (B) — P„ BPa+Pa^BPa ■ One gets that the relations 
equivalent to (|14l) . (|15[) hold with indices ^' at the relevant operators. 
The point of this construction is that it gives a 

Theorem 2.2 [Second Order Adiabatic Theorem] Under the hypotheses of Theorem \ 2.1\ we 
have for any (t, a) € [0, 1] x [0, 1], 

U a (t) = wW(t)$V(t) + 0(e 2 ), (35) 

where the error term is uniform in a £ [0, 1] . 
Remarks: 

i) If the Hamiltonian is regular enough in the t variable, it is possible to get arbitrary order 
adiabatic theorems. 

ii) These adiabatic theorems further yield the perturbative estimate 

W^it^Ht) = W a (t)Q a (t) + 0(e). (36) 

iii) This result can be found under various guises in [24], [4], [18], [25], .... As such, it is stated 
and proven in [16] , Theorem 3.3.1, p. 38 (for a fixed). Again, the uniformity in a of the error 
term is easily checked. 

As a consequence of this second Theorem, we have the approximation 

Q a (t) = - f ^ ) ~\s)Wi 1) ~\s)d a H a (s)Wi 1 \s)$i 1) (s)ds + 0(e). (37) 
6 Jo 

We point out again that the analog of the formula above with operators W a and $ a in place of 
Wi 1 ' and ^ yields an error term of order one, instead of O(e). 

We are interested in the matrix elements of Q a (t) with vectors belonging to the spectral 
subspace P a {Q)T-L, so that from now on, we focus on the operator P a (0)Q a (t)P a (0). The goal 
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is to express the information in (|37|l in terms of more familiar quantities, like dynamical phases 
and geometric phases, in certain cases. 

The next technical result says that we can express P a (0)Q a (t)P a (0) as a leading term (of 
order 1/e) that corresponds to the replacement of U a (t) by its second order adiabatic approx- 
imation with P^\o) in place of P a (0) and d a H^\t) in place of d a H a (t) plus a term of order 
e° which will give rise to the geometric contribution: 

Proposition 2.1 Under the hypotheses of Theorem \2. 1\ we have for any (t,a) £ [0, 1] x [0, 1], 
P a (0)Q a (t)P a (0) (38) 
= P a (0) (- f ^ ) ~\s)w£ ) -\s)p£\s)d a H£\s)p£\s)w£\s)^\s)ds\ P a (0) 



+P a (0) (^J\~ 1 {s)W- 1 (s)d cl K cl (s)W cl (s)$ cl (s)d s ^ P a (0) + O(e). 

Proof: Plugging the relations 

P a (t) = P a (t) 2 = P a {t){Pi 1] (t) - e]l a (K a )(t)) + 0(e 2 ) (39) 

and 

d a H a (t) = d a H^(t) + ed a K a (t) (40) 

into the expression (|37[1 . and making use of the properties of Wi 1 ^ and together with (|36[) . 
we get 

P a (0)Q a (t)P a (0) 

= P a (0) (- r^\s)wi 1) '\s)Pi 1 \s)d a H^\s)P^\s)Wi 1 \s)^\s)ds) P„(0) 



+P Q (0) (J $- 1 (s)W-\s)d a K a (s)W a (s)^ a (s)ds^ P a (0) 
-P a (0) (n a (K a )(0) J\- 1 (s)W- 1 (s)d a H a (s)W a (s)^ a (s)ds 

+ J* $-\s)W c :\s)d a H a (s)W a (s)<& a (s)dsK a (K a )(0)J P a (0) + O(e). (41) 

We want to show that the last two terms are actually of order e, by integration by parts. We 
consider the last term only, since the previous one can be dealt with in a similar fashion. Using 
property (|31[l . the integrand of this term is 

P^O)$~\s)W-\s)P a {s)d a H a {s)K{s)W a (s)$ a {s)K(0). (42) 

Differentiating the following identity with respect to a 

H a (s) = P a (s)H a (s)P a (s) + P^(s)H a (s)~K(s) y (43) 

we get 

P a (s)d a H a (s)'K(s) = P a (s)H a (s)d a P a (s)'K(s) + P a (s)d a 'K(s)H a (s)'K(s) 

= P a (s)[H a (s),d a P a (s)\K{s). (44) 
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Hence, together with (|lip . we can write 

I© = P a (0)$-\s)W-\s)[H a {s),d a P a (s)]W a {s)® a (s)K(0) 

= PaCOjO-^aJtW-^flJiTaWWaW, W' 1 {s)d a P a (s)W a (a)] $ a (s)PT(O) 

= - l£ a s (p a (o)<&- 1 (s)vy- 1 (s)a a p Q (s)w ce (s)$ a ( s )^:(o)) 

+ieP a (0)-l»- 1 (s)a 3 (^ 1 (s)a a P Q (s)W a (s)) $ a (s)7£(0). (45) 

Thus, since W / a (s) and P a (s) are independent of e, and $«(s) is unitary, when integrated 
between and t £ [0, 1], this yield a contribution of order e. I 

Remark: We can get an alternative expression for the zero'th order term by making use of the 
identity 

P a {s)d a K a {s)P a (s) = iP a {s)d a [P a (s),P a {s)]P a {s) (46) 
= iP a (s)[P a (s),d a P a (s)]P a (s). 

We want to focus now on the situation a a (t) — {E a (t)}, that is when P a (t) corresponds 
to an isolated eigenvalue E a (i), not necessarily simple, possibly associated with an infinite 
dimensional spectral subspace. That means that 

Ha(t) = £«(t)P«(t) +K(t)H a (t)K(t) (47) 

(H^-z)- 1 = P f ] + PZ{t){H a {t) - zr'Kit), (48) 
E a (t) - z 

where the reduced resolvent P a (t)(H a (t) — z)~ P a (t) is holomorphic for all z's inside the loop 
7 of the definition © if P a (t). 

This case allows to distinguish nicely contributions from the "dynamical phase" and from 
the "geometric phase" in the usual adiabatic language, in the periodic case. This comes as a 
simple 

Corollary 2.1 Assume a a (t) — {E a (t)}, then 

P a (0)Q a (t)P a (0) = P Q (0) Q J d a E a {s)ds + iW-^tfdaWcit)) P a (0) + O(e). (49) 
Remarks: 

i) We do not use periodicity in any of the variables yet. 

ii) The form of the order zero term is similar to ©, which allows to interpret it as the geometrical 
charge transported in the adiabatic process. This is supported by the fact that this term is 
determined by the spectral projector P a (t) only. 

Proof: We will approximate P& (t)Ha (t) by E a (t)Pa (t) by perturbation theory in e. By 
means of the Riesz formula © for Pj, (t) we can write 

p£\t){H£\t) - E a (t)) = -±J(Hg\t)-z)-\HP(t)-E a (t))dz 

= ~ j \z - E a (t))(H^(t) - z)- 1 dz. (50) 
We used f Idz — 0. Introducing the perturbed resolvent (|28p . we get 

P^ (*)(#«(*) - E a (t)) = -± J{z - E a (t))(H a (t) - z)- 1 dz 

~il I (Z ~ S <*(*))(#«(*) - zy L K a (t)(H a {t) ~ z)- 1 + 0(e 2 ), (51) 
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where the remainder keeps being of oder e 2 when differentiated. By making use of (|48p and the 
fact that the reduced resolvent is analytic inside 7, one gets from Cauchy formula that the first 
term of the right hand side is zero whereas the second yields 



f{z- E a (t))(H a (t) - z)~ X K a {t){H a {t) - zy 1 = eP a (t)K a (t)P a {t). (52) 
J 7 



e 

This term is zero due to (|16p . hence 

P^\t)H^(t) = E a {t)P£\t) + 0(e 2 ), (53) 
where the remainder term can be differentiated. Therefore 

d a H^(t) = d a (E a (t)Pi 1) (t) + P^(t)Hi 1) (t)'K(t)) + 0(e 2 ) (54) 

and we get 

P^(t)d a H^(t)PP(t) = P£\t)d a E a {t) + 0(e 2 ). (55) 

This allows to further simplify the first term in the expression of Proposition 12.11 making use 
of (J29J) and (31]) to get 



P«(0)(-/ \s)wi 1) \s)Pi 1 \s)d a Hi 1 \s)Pi 1 \s)Wi 1 \s)& a 1) (s)ds)P a (0) 



= p ot (Q)- f a Q s Q ( s )Pi 1) (o)dsP Q (o) + o( e ) 

6 Jo 



= P a (0)- [ d a E a (s) ds + 0(e). (56) 
e Jo 

The last term in the expression of Proposition 12 . 1 1 is dealt with as follows. The condition (I47|l 
implies 

$?(t)P«(0) = P a {0)e*i- r ° E °-( s)ds (57) 

so that we are left with 

P«(0) (J W- 1 (s)d a K a (s)W a (s)ds^j P»(0). (58) 

The argument leading form Q to ((5]) depends on the differential equation satisfied by U a (t) 
only, and thus applies to W a (t) as well, whose generator is K a (t), mutatis mutandis. This ends 
the proof of the Corollary. I 

Further specializing to the periodic case we get 

Corollary 2.2 Assume a a (t) = {E a (t)} and suppose t i-¥ H a (t) is periodic in t, of period 1. 
Then, 

i) if E a (t) is non-degenerate 

1 f 1 

(<p a (P)\Q a (l)<p a (0)) = - / d a E a {s)ds + d a p a + 0(e), (59) 



where ip a {0) is any normalized eigenvector at t — and /3 a is the corresponding geometric or 
Berry phase, 
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ii) if E a (t) is degenerate and {ipa\o),\ r £ N} denotes an orthonormal basis of P a (Q)T-L, we 
have 

<vi r) (0)|Q a (l)</£'(0)} = - [\ a E a {s)ds + i{^\t)\d a ^:\t))\l + 0{e) > (60) 
e Jo 

where ^ (t) = W a {t)^\o), r £ N, 

Hi) if E a (t) is finitely degenerate, and if {tpa (i), \ r £ {1, 2, • ■ • , A^}} denotes a C 1 , 1-periodic 
orthonormal basis of P a (t)H, we can write 

1 f 1 



(^L r) (0)|Qa(l)^ s) (0)) = - f d a E a (s)ds 
e Jo 



+i ^>* v (i)i?r(i){v# '(o)ia Q v£ 9) (o)) + Y,Br{i)d a Br{i) 

-(Vi r) (0)|a.Vi 3) (0))j +0(e), (61) 
where B a (t) solves the ODE 

B a (t) = F a (t)B a (t) , B a (0) = I, (62) 
with T a (t) defined by its matrix elements in the basis {i[>a (0)}r=i,~-,N 

(^i s) (o)|r a (t)^(o)) = -<vi s) (t)l^ r) (t)>. (63) 

Remarks: 

i) An explicit quantity for the geometric part of the charge transported is always available in 
the non-degenerate case only, see ()25|) . In the degenerate case, the geometric part is determined 
by the solution to a (second order at least) ordinary differential equation. No explicit solution is 
available in general and, moreover, the equation is parameter free which forbids an asymptotic 
analysis. However, as we explain below, there are special cases of interest in which an explicit 
expression is available for this geometric contribution. 

ii) The third point is a mere restatement of the second one, making use of an a priori time 
dependent basis of the eigenspace provided by an independent spectral analysis. 

Proof: To get the second statement, we compute 

i{^H0)\W-\t)(d a W a (t))ip^(0)} = (64) 
i{^\mW-\t)[d a W a {t)^:\Q)) - WZ\t)W a (t)d a <p$(0))) = 
i{<p£\t)\d a <pi s) (t)) - i{<pP{0)\d a <p&(0)). 

The first statement follows from (p^ (t) = ip„ (£) = <p a (t) together with the expression (|21[) . 
Finally, the third statement is proven as follows. Let us introduce 

^ r) (*) = W„(t)V£ p) (0) and ^(()=KW< r) (0), r=l, — ,N, (65) 

which defines the unitary V a (t). The link between these two bases will be made by means of 
the unitary operator B a (t) defined by 

B a {t) = V~\t)W a {t). (66) 
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By construction, [B a (t),P a (0)] = for any t £ [0,1], so that [B a (t), P a (0)] = as well. We 
compute 

B a (t) = B a (t)W-\t)W a (t) + V- 1 (t)V a (t)B a (t), (67) 
where the first term of the right hand side is zero due to (|16[) . Hence, introducing 

r o (t) = V r a _1 («)Va(t) = -V^ifiVvit) (68) 

whose matrix elements in the basis {ipa (0)} r g{i,2,---,JV} read 

(^(o^r^^o)) = -(V^(t)|^ r) (*)), (69) 

we get that B a (t) is indeed determined by the ODE 

B a (t) = T a (t)B a (t) , B a (0) = I. (70) 
Writing W a = V a B a , we compute 

(^\t)\d a ^\t)) = (^(0)\d a <pi s) (0)) (71) 
+ (4>L r) (0) IB' 1 (t) V- 1 (t) (8 a V a (t))B a (t)i>^ (0)) 
+ (V'i r) (0)| J B- 1 W(9 Q B a (t))^ ) (0)). 

With the short hand B% r (t) = (^£ ?) (0)|£ a V>£ r) (0)), we have 

S a (i)V»£ s) (0) = £>£ ?) (0)Br(i) (72) 

K,(t)B a (#i s) (0) = J^S^WW, (73) 

and 

(d a V a (t))B a (t)^(0) = ^B«' s (t)((a a T4(*)V»i 9) (0))-K(t)5 a Vi 9) (0)) 
= ^>r(t)(&^i 9) (t) - 14(t)d Q ^>(0)) 

(a a s a (t))vi s) (o) = ^d^Brm^m-B^d^io) ^ 

Inserting these expressions in (J7TJ, we get 

(^(tJISaV&H*)) - <^ r) (0)|9 Q ^ s) (0)) 



-<^ r) (0)|9aVi s) (0)>, (75) 
which yields the result. H 

Note: The operator B a (t) and its generator T a (t) depend of course on the choice of or- 
thonormal basis {ipa (t)} r =i,--,N- It is not difficult to check that if one makes another choice 
of orthonormal basis {Xa {t)}r=i,—,N such that Xa\t) = S a (t)il)a\o), then the corresponding 
generator denoted by E a (t) is related to the previous one by means of C a (t) — S~ 1 (t)V a (t) 
according to 

£«(*) = C a {t)r a (t)CZ 1 (t) + C a (t)C* 1 (t). (76) 
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3 Example 



We consider here an explicit class of Hamiltonians which, on the one hand, display permanent 
degeneracies, and, on the other hand, allow in some cases for explicit computations. Moreover, 
the physical situation considered in |10] is governed by a Hamiltonian of this class. 

Let {zi, zi, ■ ■ • , z n } be a set of n complex numbers, which we denote by the vector z — 
(zi, ■ • • , z n ) T G C™, and let £ £ I. Let us denote the standard scalar product in C n by { ■ | ■ }. 
We consider the self-adjoint Hamiltonian 

H(z) = 



( E ft 
zi 



Zn \ 









E (z\ 

\z) 



in+l 



(77) 



0/ 



relative to the canonical basis {eo, ei, • • • , e n } of C n+1 . We made z explicit in the notation 
because these parameters will become time-dependent below. 

If z 7^ 0, the rank of H(z) is equal to two, so that its kernel if of dimension n — 1, for any 
value of the parameters. If z = 0, the kernel of H (0) is of dimension n. Actually, it is easy to 
see that 

a{H{z)) = - y/E* + 4|[z|P), 0, \{E + y/E> + 4||*|[")} , (78) 

where ||z|| 2 = Y^j=i \ z j\ 2 > an d where the eigenspace corresponding to the (n— l)-fold degenerate 
eigenvalue is given by 

ker(ii"(z)) = 



( a ° \ 



We can rewrite with a — (oi, 



a n ) T e 



6 C n+1 s.t. ao = and ZjCij = 

3=1 



ker(H(z)) = © {a G C n s.t. (z\a) = 0} = ( 



(79) 



(80) 



where z denotes the orthogonal of the vector z G C™. It is now easy to express the projector 
P{z) on the degenerate spectral subspace ker(_ff(z)) in C n+1 . Let z = z/||z|| G C n and \z){z\ be 
the projector on the vector z in C™. Hence, 



\z){z\ 



is the projector on z ± in C n . Thus, expressed in block diagonal form in C n+1 
write P(z) as 

'0 

,0 P ± (z) 

Hence, with the same notations, P(z) can be written as 



P{z) = 



P(z) 



1 

\z)(z\ 



so that the range of P(z) is generated by the orthonormal basis 



Ran P(z) = span 



: span {eo, z} . 



(81) 

:©C", we can 
(82) 

(83) 
(84) 
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Let as assume now that z = z(t) is time-dependent, in such a way that [0, 1] B t z(t) G C' 1 
is C 3 . By changing the phase of z(t) if necessary, we can assume 

It is now straightforward to check that the parallel transport operator W(t) is generated by the 
self-adjoint operator K(t) = i[P(z(t)), P(z(t))], with 

K(t) = i(\k(t))(m\-\m)(m\)- m 

Recall that z = (0, z) T is a normalized a vector of C n+1 . Note that condition (|85l) is equivalent 
to saying 

z(t) = W(t)z(0) and ( J ) = ^(i) ( J ) • (87) 

Thus the determination of W restricted to P is complete. With these preliminaries behind 
us, we can turn to the interesting task from our point of view, i.e. the determination of W 
restricted to P. From (184[) above, it is clear that we can restrict attention to C" ~ eg, where 
RanP(z) Ueo ~ CI. 

Let {(fx, tpi, • ■ • , cfn-l} be an orthonormal basis in C" of z(0) ± . Then, using the same 
notation for tfj G C™ and (0, l ~pj) T G C n+1 , we have for any j = 1, • • • , n, and any t G [0, 1], 



<p j (t) = W(t)<p j €C 



l p j (t) = -\z(t))(Kt)\n(t)) eC» 



Actually, computing the parallel transport operator W(t) restricted to ker(H(z(t)) for the 
model (|77|) where z(t) G C n is given, amounts to determining n — 1 vectors (fj(t) in C™ such 
that for all j, k = {1, • • • , n - l} 2 and all t G [0, 1], 

(Vi(t)W)) = (89) 
(*>j(t)Mt))> = ^fc (90) 
(W(t)|^*(t))) = 0. (91) 

Indeed, if (|88l) is satisfied, the conditions above are met. Conversely, if the first two conditions 
above are satisfied, we get that {z(t), <pi(t), • • • , Vn-i} form an orthonormal basis, for all t's. 
Moreover, the third condition implies that tfj = co(t)z(t), for some coeffcient co(t) G C. Differ- 
entiation of (<pj(t)\z(t)) = yields ao(t) = — {z(t)\tpj(t)), so that equations (|88p are true. 

Eventhough the generator of W restricted to P(z(t)) is rather simple, these equations cannot 
be explicitely integrated in general. We present some special cases of interest which allow for 
explicit formulas. 



3.1 Special case 

We consider here a special case for n = 3 that is of interest for the physics of charge pumping, 
[10| . Let us consider the Hamiltonian 



(E 


zb 


Zl 


*\ 


zo 











Zl 











u 








0/ 
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in the canonical basis. We assume that 

\ Zl \ 2 + \z 2 \ 2 > 0, (93) 

so that a set of normalized eigenvectors corresponding to the degenerate subspace of energy zero 
is given by 

l^iy = iVi(22e 2 -zie 3 ) (94) 
IV^y = N 1 [(zf+zi)ei-zo(zie2+Z2e 3 )], (95) 

with 



JVi = i/v / kTFTNF (96) 

N 2 = l/^k? + ^ 2 2 | 2 + ko| 2 (ki| 2 + k2| 2 ). (97) 

We now compute the differentials of these eigenvectors, in order to get the generator of the 
non-abelian transformation. At this level, we allow all parameters to vary, with the condition 
that (|93[) holds. Straightforward computations yield the (negative of) the matrix elements of 
the matrix T^, , with respect to this instantaneous basis of eigenvectors of ker H 

{ipi\dipi) = iNilm (zidzi + zidzi) (98) 
{ip2\dip 2 ) = iNilm {2(z 2 + z 2 )(z 1 dz 1 + zidzi) + \z \ 2 (zidz~i + zidzi) + (\zi\ 2 + \z2\ 2 )z dz } 
(ip2\dipi) = N 1 N2Z (z 2 dz 1 - z\dz2) 
{ipi\d4> 2 ) = -{tp2\dipi). 

We can simplify some more this matrix by passing to the time-dependent basis 

Xiit) = ^( 7 (t))e l/3j(t) , Pi{t) = i ( mdipi)(l{s))ds, (99) 



where the integral is taken along a path [0, 1] 3 t >-¥ 7(t) in the parameters space. The matrix 
F x corresponding to the basis {xi, X2} °f eigenvectors of ker H now reads 



-e i(ftj - w (V>i|#2) 



(100) 



Setting 

x(t) = e- l(/32(t) -' 9lW) (V'i|dV'2)(7W), (101) 



r xW = ( JL n W ) > ( 102 ) 



so that 

r f+\ — 1 

x{t) 

we have to solve the ODE B(t) = T x (t)B(t), see (|70[) . to determine in general, no explicit 

solution to (I70|) with such a matrix can be obtained. 

However, in case x(t) = p(t)e*®, where ■& is constant in time, T x (t) = p(t)M, where M = 

, ( , and B(t) is explicitely given by 

U Vsin(/ ( >( S )d S )e rf cos(/ t p( S )d S ) ,/ V ' 

We consider below a case of this type, which allows to determine explicitely the geometric part 
of the transported charge over a period. Moreover, we express the geometric content if the 
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parallel transport within the permanently degenerate kernel of H as a solid angle in the space of 
parameters, in a similar fashion to what is done for the Berry phase, in case of non- degenerate 
eigenvalues. 

Let us assume that 

Zj — e' 9j rj and dzj = e' 6j drj, (104) 

that is, only the moduli of the complex numbers Zj vary with time. Plugging this into (|98l) 
yields 

= (105) 
fatal #2) = 2iNl sin(2(6»i - 2 ))rir 2 (ridr a - r 2 dn) 
fata|#i) = iViiV 2 r e l9 ° (e-^-^Wn -c i(9l - 92 V 2 dr 2 ) 

with 

iVi = , 1 iV 2 = 1 (106) 

VR + rl)r 2 + (rf + r\ + 2r\r\ cos(2(0! - 2 ))) 

Further assuming 

9i = 02 = 0, (107) 

we finally get 

r *- (r?+ri)Vr 2 +r 2 + r2 4 e «o o J' (1 ° 8) 

which is of the form (|102|l . The argument of the sines and cosines in (|103[) after a period 
caracterized by a loop 7 in the space of parameters is denoted by 

Q= [ -rofridrs-rsdn) 
A (rf + r|)V^ + r?+r-| 



so that 



cos(fi) sin(f2)e~ iflo 



B W = ( k - Hi a(n)> com )■ (110) 

Similarly, if Z\ and z 2 are as above and zo — tie 100 + i 2 with t%,ti real and dr± = dr 2 = 
d9o = 0, with + r\ > 0, we have 

IV =isin(fl )(ti<tta-tactt l )/(*i+*l+ 2*1*2 cos(flo)) °V (111) 

3.2 Geometric interpretation of Q 

The explicit computation of fi possesses a nice geometric interpretation, see (|119[) . as we now 
explain. 

For notational convenience, let us introduce cartesian coordinates (x,y,z) — (ro,ri,r 2 ). At 
the end of the loop 7, we have 

z(xdy-ydx) . 



(x 2 + y 2 )\/ x 2 + y 2 + z' 2 

W+y^W^WT^ ( j ' \dz) ' (113) 

Applying Stokes' Theorem, we can replace the above curvilinear integral by the flux of the curl 
of the corresponding vector through any surface £ such that 9E = 7, oriented consistently with 
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the orientation of 7. 
We compute 



curl / X = (x 2 +y 2 + z 2 )- a/2 \y\ (114) 

= 4 (H5) 



where r = ^ y j and r = x 2 + y' 2 + z 2 . Hence we can write, 

[ -^dy-ydx) = [ L dB . (116) 

J 7 (x 2 +y 2 ) v / x 2 +y 2 + z 2 J^r 3 

Consider the projection 7 of the loop 7 on the sphere S 2 described by the unit vector f = f/r 
along 7, and define E C § 2 such that 9E = 7. Now, we can choose for E the surface which 
coincides with E, and joins 7 and 7 along rays parallel to the unit vector. Since the flux of ^ 
through the latter portions of E is zero, we finally get 

j -z(xdy-ydx) = r f ^ 

J-y (x 2 +y 2 )sjx 2 + y 2 + z 2 J^r 3 

du, (118) 



is 

32 



where we used the fact that on § , da = dujr, with dto the differential of the solid angle. 
Therefore, we have obtained 

n = fi(E) (119) 

where f2(E) is the oriented solid angle described by 7 through § 2 . If 7 is oriented positively, 
Q(E) > 0, and, in any case, < |fi(£)| < 4tt. 
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